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A reliable system for performing a Fourier analysis of a 
reactor noise signal has been constructed and utilized. Both 
auto- and cross-correlation techniques were used to study 
noise signals from the MITR hohlraum and the exponential 
lattice. As predicted by theory, the efficiency of the detec- 
tors used was insufficient to provide positive results. The 
lattice transfer function was electronically simulated and 
analyzed by the system confirming the reliability of the 
equipment . 

Several methods are suggested for increasing detector 
efficiencies to provide signal levels which will permit 
extraction of information about the kinetic behavior of the 
test lattice * 
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INTRODUCTION 

1 ol Background 

In studying the neutron characteristics of any nuclear 

reactor core assembly It is desirable that the research be 

done as economically as possible , and preferably that 

several independent methods be used to determine the 

desired information. he MIT Lattice Research Project has 

for some time studied anous reactor parameters in a heavy 

water=moderated exponential facility located adjacent to 

( ) 

the MIT Reactor o 1 ' ^ 1 See Figure 1 ol for a cut=>away view 

of the MIT Lattice Research Facility - 

In recent years „ considerable use has been made of 
the pulsed neutron technique in determining k for test 

UO 

lattices in order to provide data supplementing the more 
common static parameter measurements It is the purpose 
of this thesis to explore the utility of the noise analysis 
technique as a means of complementing the pulsed neutron 
and other methods of measuring the various parameters 
measured in the MIT heavy water lattices- 

1 -2 Noise 

The term “noise® 1 refers to th° random or stochastic 
fluctuations in the output signal from a system- Reactor 
noise is thus the random fluctuation of power level due 
to the discrete nature of the multiplication process- 
These fluctuations are generally small and iruaXly only 
about one percent of the total power - In order to 

observe this noise, a detector must be used and often 



2 



DENSE 

CONCRETE 




FUEL 
ELEMENTS 

EXPERIMENT 
PORT 



Port 3CH2 



_thermal Column 
Door 

LATTICE 

facility hohlraum 



OBSERVATION 
WINDOW 



MEDICAL THERAPY ROOM 



ATER SHUTTER 



Fig. 1.1 Cut-away View of MIT Lattice Research Facility 
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the term reactor noise Is used to describe the fluctuations 
of the output signal of the detector » It is important to 
remember, however, that the detector only seej/s a token 
sample of these fluctuations a. id that, owing to the finite 
size of the detector, the spatial-dependence of the noise 
measurement is blurred . Only recently has the space -depen- 
dence of reactor noia been closely examined,^^ Throughout 
this paper, howevei , sp ,• ^-independence will be assumed 
since this will piove au^quate lor ^je&' n .irpofc =s „ 



1 »3 Summary of Theoretical Curmid’: rations 

Since a reactor depends upon neutron multiplication for 



its operation, it is reasonable to expect that among the 
various neutrons which aie registeied by a detector, many which 
have a common ancestry wil.1 produce pulses superimposed upon 
the doC, level of the signals The number of such related 
pulses will depend upon the detector til l rry as will be 
seen later . While the detector d„c level is Itself main- 
tained by the accumulation of many pulses, we die unable to 
extract any information involving the discreteness of the 
multiplication process from it, and it thus holds no interest 
for us here , Thus one can in • hi focus attention upon a 
fluctuation In the signal due tu the daughter neutron of a 
given fission, followed a certain A- later by another 
impulse due to a grandaughter neutron of the same fission. 



Since this time interval At is a ium 
lifetime in a given system, it . - i 
information about the multiplying gi * 



1 n of the neutron 
t- mss! 1 ' to extract 
,5 ot tha medium 



by submitting the a c (or »» l ■-*) port l>'*n of the detector 
signal to analysis o One mat hod of analysis is to determine 
the autocorrelation functlortv wM.< h i-. defined as the expec^ 
tation value of the j » i d u if t,K: sirn-ii :*t a time t 9 , and 
the signal at a time r.,s 

* j 

0 nn ( t ) = E|n(t ,)n(t j\ {1.1) 

where for simpll< it Mo il i, twicer to he the neutron 

density n(t). and 7 * t F 11 im th* development of 

Bierman, et one* 1 an e -i lu-xte )^ rH ( t) by first con- 

sidering the standard reactor kint.ti 6 equations in their 
space=>inde pendent , linear, t i.we invariant f ,iras 



dn 

dt 




u 



4 



i 



X « * S 



P-2) 



dC 



1 = - A .0 . + k -- 



~ar 88 “^i w i T 0 (I ° 3 ^ 

Since only pt"nu«f n?utx *n h ha 1 * Is- Ff importance , 
these equations may he pit in the forms 



if + an = k | n o * S “ 1 ^ t 



(l„4. 



where a 






the prompt neutron decay constant, and 



the Inhomogeneous terms are cor u, f - 'ed as the forcing func= 
tion f (t) . 

The solution to Eq* (lo4) may be written as a convolu« 

tion integrals 



n 



oo 



f 00 

W(t-7.)f(X)dT. '•'ft) { • C)- (1*5) 
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where W(t“£) la tiie impulse response function of the system. 
This function Is obtained by normalizing the homogeneous 
solution to Eq « (1.4) and is found to beg 



W(t-'C) = 



o axp[- c<(t~tj| t>Z 



0 t < t 

, 

One may evaluate Eq, ( 1 1) ;.1> Ingg 



( 1 . 6 ) 



f* nn (t) = | e * at! • (1-7) 

If one were to assume that the forcing function were 
instead a burst of neutrons, as in a pulsed neutron experi“ 
ment, f(t) could be approximated by a Dirac delta 6(t) giving 
as a solution to Eq, (l,4)s 

n(t) = oe“ at , 



At this point, the complimentary nature of the pulsed 
neutron and noise analysis experiments becomes evident, for 
each provides an independent method for determining the 
prompt neutron decay constant a. 

Various noise analysis methods are available for deter~ 
mining a by studying the a .c » detector signal in the time 
domain, however, by transforming Eq« (1,5) into the fre- 
quency domain by the use of Fi’ ier transforms, a method of 
analysis which involves less data processing can be employed. 
The Fourier transform of the autocorrelation function 
is called the power-spectral“density function and is given 
byg 



V'n>> 



1 ^ 

TT 



foo 



=>00 



nn 




( 1 . 8 ) 
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which when combined v, ltn x 



tv U) = — 

t nn* 1 7 r 



cT 



F i.o W 



a 

j 



■ y; 



The spectral -fit- ■. l <■ j !*■ 1 1 ■*' a lew pass 

filter with a breakp m+ <> : 1 , 1 ! J 1 arier analy- 
sis of the dete 'tm uigua • should be possi- 
ble to detex'mine •> tne .■■ ■ ' 1 ■-> 1 tart 0 

1 0 4 Detector £i l 

Xn oi'dci fcv <i ■ * 1 : n ' non 

and noise measure me* r- , 1 ■ -■ 1 . i that- both 

experiments are hampered cv * 1 ' 'i v u,. ignal levels . In 

the pulsed neutron expei’c n*. . ■ 1 “ form of residual 

neutron levels due t < n- ut 11 in the ^vstem, and 

can be overcome by in • ' "<■* intensity 0 In noise 

experiments, as we shall . m - ti n, the detection 

process itself previd B a Pi4 ' spectrum, 

which can only be ovei-< *-**i*<- '*'* 1 1 n ' f ■ ■ ~ r el iiciency 

or by performing er-a;- c r : "-m • m These limita- 

tions can be described it.’.'n Ml 1 ‘-round as an 
additive constant t r 1 * 



N(t) = A p + B p : 



at 



V- 

+ l’° ' 



( 1 . 10 ) 

( 1 . 11 ) 



n » { ) = A + “ 

^ nr> 1 n i 

where now a figure of merit fvc m-t -/Stem can be defined as; 



JB_ 

A' + 'B 
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The utility and feaslbilit of noise experiments thus depends 
upon the magnitude of the t Igure of rre rit 

By nursuinR the analvtical method of Gohn^'^ » and extend- 
ing it to the case of a sub^riticil reactor, the effects of 
detector efficiency can be studied. The basis of this method 
rests on what is known as the S-hottky formula, which was 
originally developed t> iescribe the noise in a temperature 
limited electronic ii' t tss 

<IXI 2 > = 2e 2 m , (1.12) 

p 

where <Jlj > is the spectral density cf 'he diode current 
2 

noise in amps sec, e is the charge carried by each electron 
in coulombs, and m is the av- rage number of electrons flowing 
per second. Cohn adapts r i* " r-mula to the reactor case 

giving % 

<!S 0 I 2 > = 2 1 q 2 m J (lo!3) 



where CIS l 2 > is the spectral iensity of the noise equiva- 
o 

p ^ 1 

lent source in neutrons se , q, is the net number oi neu- 
trons produced in a nuclear process of type i, and is the 
average number of nuclear processes of type i occurring per 
second in the reactor. As ah< In Appendix A this formula- 
tion results in the following expression for the noise equlva* 
lent source in a sub critical assemc . 



<1S 0 I 2 > - J- 



2 



k - -■ + “{ V - 0 V + 1 ; 
v It 



(i.i4) 

This source is transformed by re a + ■" •<* ir - transfer 

function T {to) into the neutron s. • 



sitys 
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< j n(cJ>) l 2 > = \T(*«.>) I (\ h J > (1. oX5) 

By considering the power sne tral density of the detec- 
tor current due to b th tt > i t n un<_ or related neu- 

trons, it can be shown i- t- u- i u>-Um tor current auto- 
power-spectral-denslt y 1* 



2 , 



<|I a (^)/ > = ^ F oj 



I ^,1 

i'i ► - i -■ 



. b i 

* L 






P 



=P<* 






1 



-*T f 



Uur+as 
( 1 - 16 ) 



where 

G 



F 

k, 



s detector efficiency 



= fission r-- 1 e 

= k(i-p) 



k .n 




fissions detected 
fissions occurring 



= average number *f neutrons emitted per fission 
q = charge collected per fi - ’ i r n defected 

!tM I 2 = 

Thus * for an autot oriela *n •'tea -ur- t , the detector 
efficiency must be such that the right hand term is greater 
than unity, or for <+>« as 




Using the accepted value 






( 1 - 17 ) 



v + 1) 

q ) of 2.5 for the factor {J and 



7„47 for 7% this t. 







(1 ol8) 
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Figure lo2 displays thioi f fI *. - 0 90 o 

Since detector effic ■ , ■ - v ve 10 are v rv difficult 
to obtain without aigm.flr mtlv . 4 r urbl-.c. the system under 

study, a major stn . 1 i ^ ■ 1 11 -t noise analysis 

in far-subcritioa 1 ayst: ?unr (■ ■ v ■ . bei | 'yri- i obvious# Some of 

this problem could be all- d if the up> or related noise is 

filtered out „ By ua ' 1 * tors and multiplying the a»c» 

components of the u r 1 wfalch did not coin— 

cide would tend to c *n« , i -f •.■•*- ■ which did coincide 

would reinforce ea^b - 1 „ , - r i rf time, the 

correlated portion of the -d*ec t < uro would tend to be emphasized o 
This method of cross cor re i at , w«s n p d by Seifritz* 

Stegemann, and Vath In m -king * 'ax urn n n the coupled fast- 
thermal reactor STARK J 11 *’ rte : -e?ult . indicated that with 

their equipment, a ratio f -^ed to unc or related noise 

of Ocl was suffi cient to >1' " 'f ' * if cross- 

correlation techniques wcx* 1 ‘ 1 » bis i < a ia-'" t 'Gr of twenty 
improvement over the rati 1 ■- needed in auto- 
correlation experiments, Of t' e relative advantage of 

cross-correlation ovei au -i cm iation experiments depends 
upon the individual rd h w the signal is 

analyzed „ Figure 1 L -hr 1 • ■ -*u -related detector 

efficiencies required has- -V i-rieoce. Notice that 

detector eff iciem ies of • 1* ^ - r*' cent axe still 

required for a k of 0 ' ’ n * e ™nge of P ri ~ 

mary interest in MIT labti.-. 

1 °5 Purpose of K?es ^nt Jf ' 

This thesis ir •- < n . • • ■ : r a * *- <nt cing a 

reliable analysis system t ■ - : 11 - r aut«. 



Detector Efficiency ( 
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cross-cori'elatir? Prora the start, it 

was realized that on n <? -<u h .1 tv - vj t onitri j, the 
problem of att-ainip, suf • - '• “tor efficiencies 

would have to fc; cn - 

Detector efficient.' in - ,1 * matte’ of detector size 

and neutron sensitivity 1 ■ ■ - readily available, 

the Westlnghouse I „ f. 1 - < r hasher provided the 

best combination c£ n . • * y / ■ v h )4 amperes/ 

neutron/cm /sec) an*3 ^i*. ( t -lii-ve area), -As 

such, this detector #, . ir.ienr enough (in the 

hohlraum and lattice, < ' % i • ) ,0 • ^sare the spectral 

density; however it ^ tv >e Ri • it provided 

a noise signal compara it 5. 1 ■<=? t -> those which would be 
studied in the lattice iat: > < 1 tv us provided us with infor- 
mation on the design rate ' i - < the analysis system. 

Second, any very' la - e «i the 4i*ihlraum which might 

be present due to spaiial vlfec >e> uii r still be detected. 
Thus, both auto- and l roo- . v - ■ • ^ r . ; mentis were made 

in the lattice to insure that a reliable system would be 
available in the future, concurrently, a separate study was 
being made of methods f * 1 inv‘ 1 ■ t, d- : terror efficiencies 
to the levels needed 

Chapter 2 descrit ; t r : • ■ ■ r he initial phase 

■ aeirJbly aiv ... s> r t ¥ „ Chapters 3 

■ni . . ./ both 

-iv .** ( pter 5 



i 
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CHAPTER 2 

ASSEMBLY OP EQUIPMENT AND PRELIMINARY 
HOHLRAUM MEASUREMENTS 

2,1 Assembly of Equipment 

Several methods of analyzing noise signals have been 
developed, often classified as to whether the analysis is 
done concurrently with reactor operation, or done afterward 
by analyzing either an electronically or mechanically (pen- 
recorder) recorded signal, Whether the analysis is done 
digitally or continuously, further specifies the method of 
analysis to be used. The choice of method of course depends 
upon financial and/or time limitations, and also upon the 
degree of accuracy desired. 

For present purposes. It was decided that the use of 
analog computing equipment would provide the most flexible 
system of analysis. Not only would it provide a real-time 
method of continuous analysis but it would also be adaptable 
to other uses such as reactor dynamics simulation. For this 
reason one George A, Philbrick MU/DV Duplex Multiplier/Divider 
and three G.A.P. Model HK Operational Manifolds (each con- 
taining five chopper-stabilized operational amplifiers), were 
acquired. The latter afford the capability of performing a 
total of fifteen separate linear operations on one or more 
signals, whereas the former performs the non-linear operations 
of multiplying or dividing. (See Appendix D for a brief 
discussion of the characteristics of operational amplifiers 
and a description of the basic cir^ used In. performing 

the linear operations.) The method of a -i sis used was to 
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take the Fourier integral < t . < , n function 

with zero delay time,, A felnok dUi. ram <i the ah trliry is 

shovm in Figure 2 . 1 . 

After using the heter. aynlnft baoht j ' multiplying 

the detector noise by a sine wave of known frequency, a low 
frequency band-pass filter of fix -q bain is used to 

pass only those noise component! differs from 

that of the oscillator i y i i* than the fix tor bandwidth. 

The resulting signal la » •• .j - , . , va ion fc, 

the integration time t , . i , , \ •pac.tral 

density at a frequency oj ( ~ in a [width A co 0 

Scanning the spectrum is all mf utei 1-v \>ir i- g the oscilla- 
tor frequency oj Q over the v 

2 « 2 Preliminary Experime nt 1 r 

Initial assembly of tf ■ ,1 i* * . . i ■ ■ i , tat ion 

were done by the author tog 

As stated earlier, this firs* pha* i e two purposes. 

First to determine if the as 1 m . i *' s p i forming 

correctly, and secondly to determine if ft- thermal column 

and hohlraum served as a true white noise source of neutrons. 

This latter question is imp.. * - jine- , » n b • seen by 

Eq. (1.15), determination nf lattice h'teter function 

requires foreknowledge of th* >p •: 1 1< • «r f the source. 

Only when the external source is itself random may the 

f i iq ) 

simplified pile noi3e formulation ! 

For these measurements, W* ivi 6377 ioniza- 
tion chamber, was place! i,i t * r-igh a hole bored 

through the graphite wall 1> 
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Fig o 2 d Block Diagram of Auto»Spectral Density Analysis System 
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Preliminary measurements showed that the total d„c. 
detector output was 175 M-amps with a noise component of 
about 1 pamp. Since it was desirable to have a voltage 
input to the vacuum-tube equipment on the order of one-half 
volt or larger, the need for a high gain preamplifier became 
apparent . To avoid excessive 60 cycle pickup and to intro- 
duce as little instrument noise as possible, it was decided 
to use an Analog Device Model 106 Solid State A o 0. opera- 
tional amplifier, powered by four F4BP 6 volt batteries „ In 
addition to meeting technical requirements, the amplifier 
has the added advantage of economy costing $ 26 „ For similar 
reasons, the detector high voltage was supplied by 10 67 o 5 V 
batteries. To avoid saturating the preamplifier with the 
do Co level of the incoming signal while still obtaining maxi- 
mum amplification of the noise component, a buckling voltage 
controlled by a potentiometer was inserted just before the 
input o Also, since there was a large low frequency variation 
in the neutron flux, apparently due to control rod motion in 
the reactor, a high pass network (Rg and Cg) with a break- 
point at 2.5 cps was introduced (see Figure 2 0 2) 0 Thus the 
preamplifier response iss 



e 



o 



-21^(1 + Tffig C 2 8 ) 

<2*2 ^H 1 + 



- i 



) 



in 



( 2 . 1 ) 



which for large s = j becomes 



R, 



in 



( 2 . 2 ) 



Although a high frequency roll-off due to finite space charge 



16 



R 1 = 510 Kfl 
R 2 = 12 Kil 

C 2 « 220 ixf 
R 3 = 100 KSL 
P 1 = 500 KXl 
E = ->12 V 
A = 150,000 




Plgc 2.2 High Gain A.C. Preamplifier Circuit 
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collection times was observed, the preamp-detector frequency 
response was constant up to 5*000 cps for a detector voltage 
of 720 V, 

The noise amplitude at the preamp output was about one- 
half volt, which was sufficiently large compared to the 
0o 01 V self-noise of the preamp „ Two stage amplification with 
a total gain of 200 provided a 100 V peak-to-peak signal at 
the input to the heterodyning multiplier. Here, the signal 
is multiplied by a 100 V peak-to-peak sine wave of known 
frequency. Since the output of the multiplier is scaled by 
100 V, the resultant signal is 100 V p-p. The combination of 
the two filters produced a 0,88 cps pass-band centered at 
0,128 cps. The inclusion of an amplifier with a gain of 20 
between filters maintained the signal at a sufficiently high 
level to make the effects of component drift at this stage 
negligible. After squaring, the filtered signal was integra- 
ted using an averaging or "leaky 11 integrator, with a time 
constant of 100 seconds. The circuit at this point is shown 
in Figure 2,3, Results of the preliminary measurements in 
the hohlraum are presented in Figure 2,4, Two features are 
important. The first, already described, is the high 
frequency roll-off as a function of detector voltage. The 
second is the low frequency erratic peaking in amplitude. 

To investigate the latter, the effect of the reactor 
source neutrons was considered. Since these neutrons orig- 
inate in the critical MIT Reactor, a low frequency breakpoint 
in the reactor neutron spectral dens*. would be expected at 
^ = ( 3 /J, , This value is approximately ,0075/, 0015 for the 



<?»OsC 
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Fig. 2.3 Circuit Diagram for Auto-Spectral Density Measurements 
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MIT Reactor * resulting in a break at f = * /2w ~ 0,796 cps. 

In order to determine the amount of correlation intro- 
duced into the lattice by these source neutrons* one can 
assume that the lattice is a detector and apply Eq, (l„l 6) 0 
To do this* the efficiency of the lattice for reactor neu- 
trons must be estimated. Less than 0,02$ of the reactor 
neutrons pass the hohlraum. This value may be taken as 

an upper limit for the efficiency of the lattice. Using 
Eq, (l„l 6 ) for the case of a critical reactor (ki^ l)* with 
a = 4,8 sec * it is seen that above about 5 cps* correlation 
of neutrons in the lattice should be negligible. Since the 
frequency range of interest in the lattice is about 3*000 
cycles* the effect of the low frequency correlation may be 
ignored. Thus the hohlraum may be considered to be a white 

noise source of neutrons for the lattice, 

-tV->e 

For this reason* t b4* behavior of the spectral density 
curve of Figure 2,4 below 100 cps caused some concern, 

2,3 First Phase Conclusions 

In order to check on the reliability of the system* 
measurements of the spectral density were taken several times 
on different days. Above 50 cps* the results were found to 
be reproduceable within the predicted experimental error (see 
Appendix D)„ Scattering of data points below 50 cps cast 
doubt upon the reliability of the system at low frequencies. 
Preliminary investigations were concluded at this point 
with the next objectives being tot 

a) attexnpt to improve detector efiu i^ncy 
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b) investigate more thoroughly the low frequency 
. response of the system 

c) investigate frequency response in the lattice 
itself. 

The problems of developing a detector with increased effici- 
ency for lattice noise analysis is the subject of another SM 
thesis by M.G. Johnson, currently underway. The next 
chapter describes work done by the present author on the 
items (b) and (c). 






CHAPTER 3 
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LOW FREQUENCY TESTS AND AUTOCORRELATION MEASUREMENTS 

Section 2.3 described the erratic behavior of the analyzer 
at low frequencies. It was necessary to determine the cause of 
this behavior before proceeding to lattice autocorrelation 
measurements. 

3. 1 Introduction 

Up to this point, the oscillator signal had been provided 
by first a General Radio Unit R-C Oscillator Type'1210-C and 
later a Hewlett Packard Test Oscillator, Model 65OA, the latter 
providing the signals of lowest frequency down to 10 cps. In 
order to look more closely at the low frequency response, a 
Waveforms Model 403B oscillator was used, providing signals as 
low as 1 cps. The filters were also modified as shown in 
Figure 3.1. 

This time the spectral density of the hohlraum noise was 
investigated in the range 1 to 200 cps, with the results shown 
in Figure 3.2. As can be seen, results were again erratic, 
with a general trend of higher amplitudes at lower frequencies. 

The erratic low frequency response suggested the possibil- 
ity of component drift, while the Increase in amplitude at 
low frequencies suggested that the d.c. portion of the detec- 
tor signal was not being completely cancelled, allowing the 
carrier wave to be passed to the filters. Referring to the 
latter problem, note that a one-half volt d.c. signal at the 
noise Input to the multiplier, when multiplied by a 100 V 
peak-to-peak carrier signal, of frequen, will result in 
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a one-half volt peak-to-peak sine wave at the multiplier 
output, At high frequencies (such that oi is much greater 
than the center frequency of the filter) this component is 
filtered out leaving only the signal due to the a,c, (noise) 
portion of the signal. However , since the filters are not 
square (effectively 2k db /octave attenuation as seen at the 
averaging circuit input), at low frequencies this sine wave 
would not be attenuated enough to be negligible. Total cancel- 
lation of the d,c, level of the signal was insured by adding 
a blocking capacitor in series with the first amplifier. 

Multiplier drift (rated at +1$) which was measured to 
be as much as 2 ¥ in one hour, accounted for the erratic 
behavior at low frequencies. This drift was minimised by 
adding a ventilating fan to provide better temperature stabil- 
ity! however, drifts of up to 0,75 V per hour were still 
encountered. Figure 3«3 shows how these measures served to 
reduce these problems, 

3,2 Simulation of Lattice Measurement 

To Insure that the equipment was indeed operating 
properly, it was decided to simulate a lattice measurement 
by adding a capacitor to the preamplifier, in parallel with 
the feedback resistor. Neglecting the low frequency filter 
network, the preamplifier has a gains 

R*. 

_ . . f. (3 i) 

i ” 1 + R^G-s 0 } 

n ll 

With R r = 51 OK (1$) and C f = . >0'. pf , (ma-asurea) this results 
in a simulated breaKpoint at; 
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At the same time* the filters were modified to widen the 
pass band* which would allow shorter integration times (see 
Appendix P). The filter modifications are shown in Fig. 3.4. 

Results of the simulated lattice measurement are shown 
in Flgo 3o 5o The breakpoint occurs at that frequency corres- 
ponding to a 6 db decrease in the squared signal and is 
estimated to be 320 + 20 cps* which is in good agreement with 
the predicted value of 314 t 4 *:p „ 



3c 3 Lattice Autocorrelati on Measurements 

Having been satisfied that the system was now functioning 

properly, attention was turned to making a test run in a 

B o 0-moderated* uranium- fueled lattice. Estimates of k ^ . and 
2 eff 

a for the lattice installed in the exponential facility during 
this thesis work are developed in Appendix C* indicating that 
a detector efficiency of around 30$ would be required to 
enable determination of a for this particular lattice,. None- 
theless* since it would be several months at least until a 
more reactive lattice was scheduled for installation* it was 
decided to make a test run to insuxv that the equipment would 
be ready for use. These measurements were made by suspending 
the 6377 boronlined detector in a lead-weighted aluminum tube 
in the center of the lattice core. In order to do this, the 
seven center fuel rods had to be removed from the 211 rod 
lattice. Since the multiplication factor of the assembly was 
already so low as to preclude measur of kinetic effects* 

the small additional decrease in react tv s, caused by removing 
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Figo 3 <>5 Experimental Determination of Simulated Breakpoint 



30 

the rods was not considered significant „ First attempts at 
taking measurements revealed that the 75 foot cable from the 
detector to the preamp caused excessive oscillation in the 
preamplifier , Although this oscillation could be eliminated 
by use of an RC compensating circuit * the signal strength was 
then so small as to be comparable to the preamp self-noise. 
Therefore* the preamplifier was placed on the lattice top* 
eliminating the need for both the long input cables and the 
compensating network. Such an arrangement produced lattice 
noise signals five times greater than the preamp self noise. 
With the equipment set up in this manner* the results shown 
in Figure 3.6 were obtained. Again* the white-noise spectrum 

i 

due to the low detector efficiency is evident. 

No further work was performed on autocorrelation measure- 
ments because the cross-correlation technique to be described 
in the next chapter appeared to hold more promise. 
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4,1 Theory 

As mentioned in Chapter 1*. cross-correlating the 
signals should decrease the detector efficiencies needed 
for noise measurements. The crosscorrelation function 
is defined as the expectation value of two independent 
signals % 



E 




(4*1) 



Although 0 ^ and 0 2 consist of noise signals super= 
imposed on doC, levels, we need consider only the noises 



^(t) - U x (t) + C 1 (t) (4.2) 

J* 2 (t) = 0 2 (t) + C 2 (t) » (4.3) 

where U and C are uncorrelated and correlated noise, 
respectively. 

The crosscorrelation function of these two signals 
is given bys 

Epi(t¥ 2 (t)] - 

+ E[u 2 {t;)C 1 (t;)] + 

(4.4) 

Since the expectation value of the multiple of two 
uncorrelated functions is zero, all uncorrelated terms 
disappear and? 

E gl(t>)* 2 <*>] = E^tJC^t)] . C^«5) 

Thus, corresponding to Eq, 1,16 for autocorrelation 
experiments, one can write for cross-correlation experiments? 






<II ( 



>2 - 2 W 0 Q 2 







p ZF +a 



where now £ ± and refer to the efficiencies of the 
separate detectors. This would seem to imply that measure^ 
ments could be made Independent of detector efficiency* as 
long as the analyzer output signal is sufficiently large . 
However* as stated earlier* Se i/fritz, et al„ found that 
cross-correlation decreased their efficiency needs only by a 
factor of 20 . This is due to the fact that the cross-terms 
in Eq. (4.4) are not totally uncorrelated . A factor of 20 
reduction in efficiency requirements* however* is a highly 
desirable* if not essential* improvement if one is to make 
noise measurements in subcritlcal systems. 




The cross-correlation measurements required the addit- 
ional equipment shown within the dotted lines in Pig. 4.1* 
consisting of one additional multiplier (G.A.P. + 25 V MU/DV 
Duplex Multiplier/Divider)* 6 operational amplifiers* 
another preamplifier* and an additional Westinghouse Type 
6377 boron-lined ionization chamber. 

The two detectors were suspended in the high-flux 
region of the lattice described earlier, side-by-side* in 
separate lead-weighted aluminum tubes. Again the preampli- 
fiers were placed on the lattice top* as close to the detec- 
tors as possible to avoid preamp oscillations. Both detec- 
tors were operated at 670 V. 
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Fig, 4,1 Circuit Diagram for Cross- Correlation Measurements, See App, 
for letter designations,* All filters same. 




4 <>3 Result^’ 
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Recall that for autocorrelation measurements the filtered 



signal is squared and then averaged » The squaring process 
naturally results in a large non^z&ezr signal * whose average 
will also then be non-gsJHfJo The cross-correlation measurement * 
in that it eliminates the uncorrelated noise signals* should 
produce a signal whose magnitude relative to the auto-c or rela- 
tion signal is given by the ratio of Eq» (4.6) to Eq 0 (1.16). 
For very low efficiencies* this ratio is merely the right-hand 
term of Eq« (l„l 6 ) and iss 



£ k p 



2 



|2 

<Jl a F> (vr c<r + a' 

For 6 ^ 10“?* anc i for the lattice under study g 



. \ d > (v) z Zr 



( 4 * 6 ) 




3 x 10 



Thus since the signal level in the auto-correlation measure- 
ments was Of l 5 >- 100 V* the cross-correlation output signal 
must be < i * 2 > - 3 x lO" 2 *' V 



This signal is so low as to be 
overridden by equipment drift „ While observing this drop in 
net output confirms predictions as to the effect of converting 
from auto- to cross-correlation* the low signal made it 
impossible to extract a spectrum <, Even with cross-correlation 
techniques* therefore* higher lattice multiplications and 
detector efficiencies will be required «, 

Experimentation was concluded at this point <> Conclu- 



sions and recommendations follow in .Civ ^ * 0^5 x 9 5 ° 



CHAPTER 5 

SUMMARY* CONCLUSIONS * AND RECOMMENDATIONS 
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5 .1 Sutonafry and’ Conclusions' 

A reliable system for making auto- and cross-power noise 
spectral density measurements on a subcritical nuclear 
assembly was constructed . The system was demonstrated 
capable of providing signal analysis over a frequency range 
of 1 to 5*000 cpso 

The assembled system was tested in several applications § 

a) a signal simulating a reactor noise spectrum 
was analyzed to demonstrate successful system 
performances under optimum operating condi- 
tions (Chapter 3* Section 2) » The measured 
breakpoint frequency 320 +20 ops agreed quite 
well with the known input of 3.14 + 4 eps° 

b) The MIT exponential facility hohlraum was 
studied using the auto-correlation technique » 

None of the results (Chapter 3* Section 1) gave 
reason to doubt the assumption that the hohlraum 
provides a white source of neutrons to the test 
lattices o 

c) Both auto- and cross-correlation measurements 
were made in a lattice „ In both cases it was 
shown that lattice multiplication and detec- 
tion efficiencies must be increased before 
noise analysis can be n«ed ' extract neutron 
decay constants (Chapters 3 ar I 4)» 
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Since detector efficiency requirements vary inversely 

as the snuare of the multiplication factor of the system, 

future studies should concentrate on lattices with k ~~ as 

ef f 

close to the present lattice safeguard limit of k eff = 0„9 
as possible . Since cross-correlation of two signals 
significantly reduces the efficiency requirements, this 
method provides the greatest promise for successful analysis 
at reasonable efficiency levels » 

3 .2 Recommendations 

In order to attain the required detector efficiencies, 
the utility of gamma detectors in the lattice annulus 
should be fully explored. Such detectors would not disturb 
the multiplying system, but would detect both prompt fission 
product gammas, which escaped from the lattice core, and 
those gammas produced by neutron capture in the cadmium 
sheath surrounding the lattice „ Alternatively it may be 
possible to develop a detection system which detects the 
Cerenkov radiation produced In the moderating DgO. Further, 
the 2$ enriched UOg lattices, scheduled for study beginning 
in February 1967* should provide the capability of attain- 
ing k eff 2^0.9, which would also reduce detector efficiency 
requirements considerably. 

Finally, in future work, where large volume detectors 
must be used, the space dependence of reactor noise measure- 
ments should be carefully considered. 
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APPENDIX A 

NOISE-EQUIVALENT SOURCE FOR A SUBCRITICAL ASSEMBLY 



Neutron processes in a multiplying system are analogous 
to the random flow of electrons in a diode „ Using this 
similarity* Cohn has adapted the diode noise equation of 
Schottky to the noise in a critical reactor By S i m pi e 
modifications it is possible to extend this treatment to the 
case of a subcritical system „ 

Cohn® s adaptation of the Schottky formula iss 

<|S 0 i 2 > = 22 q i 2 m 1 , (Aol) 

where s 4IS 0 J > is the noise -equivalent source; 

q^ = net number of neutrons produced in the nuclear 
process of type 1; and 

= average number of processes of type i occurring 
per second in the system. 

Table Aol lists the processes* their rate of occurrence, 
and the net number of neutrons produced in that process o 



TABLE Aol 

CONTRIBUTIONS TO PILE NOISE SOURCE 

Average Rate Net Number of 

Nature of Process of Occurrence Neutrons Produced 



Nonproductive leak- ~ A ~-^ „ =1 

age or absorption ■*- A + 11 

Fission giving rise n __F p 

to N prompt neutrons A + F *N 



External source provid- 
ing neutrons to main- — (1 - k ) 
tain the pile -C 



1 



In this formulation delayed neutron effects have been neglee- 
tedo Here 

A = macroscopic cross-section for all non-productive 
processes Including leakage, 

P - macroscopic cross-section for fission, 
n = total number of neutrons in the system, 

A ~ Prompt neutron lifetime in sec, and 
=_//k = neutron generation time in sec „ 

A and P are subject to the neutron balance condition 

A + F p = k p ° 2 ^ 

where v is the average number of neutrons produced per 
fission « P n is the probability that N prompt neutrons will 
be produced in any one fission, subject to the conditions 
oo 

£ P»t — 1 (Ao3) 

N=1 N 



S N P = (1= 
N-l w 



V "* v 



00 (p 

2 N“P = (1* 
N=1 W 



~ V 1 " 



(A.5) 



Substitution into Equation (Alol) yields 



< is„r> = 



2n 



^AWf 



GO 

A+kF 2 | 
N=1 



■1) P N + (1= 



Applying conditions (A .3), (A«4) and (A .5) yields 



<ls o ,2 > -J- 



1 - •=• 
L v. 



4-{^-2P+l) 

y 



l-k 



(Ao?) 



The first term in square brackets is ti ^erm due to 



4o 



nonproductive (e»g., (n,/) capture and leakage) processes, 
the second is the term due to fission processes, the third 
is the term due to the external source „ Note that the exter- 
nal source term has been assumed to be frequency Independent, 
and that the entire noise equivalent source is thus a white 
noise source o Since the fission rate iss 
nk 






p ° = 57 ' 

Eq 0 (Ao?) becomes 



(A.8) 




k D “ 

+ 



2V + 1. 



(A.9) 



The source neutrons are then correlated by the multiplying 
process » The correlated component of the noise spectrum is 
obtained by multiplying equation (A»9) by the nuclear system 
transfer functions 



<Jn(f)| 2 > = I T(f ) \ <IS 0 | 2 > 



2DF 



k, 



P 






2 



■2v + y ) 



UJ 2 + 



l=k. 



(A .10) 



where a = - * is the prompt neutron decay constant , 

The average current I in a detector due to all fissions 
detected will be s 



nk. 



1 * €ci T£ =6Q3? ° 



(A oil) 



where Q is the charge collected per fission detected „ Super- 
imposed on this d.c. level are noise components due to the 
correlated neutrons and the randomne^ f the collection 
process . The former is given by Eqs » (AolG) and (A oil) 












i 

” f 



4l 




(A .12) 



(A ol3) 

The latter can be obtained by using the Schottky formulas 



For £ < < 1* these two terms are first order uncorrelated 
and are added in quadrature to give 



This result enables estimation of detector efficiencies 
needed to extract information from noise signals as shown 
in Chapter 1, Section 4. 

Note that throughout , it has been assumed that the 
analysis equipment did not introduce any spurious frequency 
dependence o Such an assumption is justified by the "white" 
noise spectrum displayed in Figure 3 <,3° 



</I c | 2 > = 2Q 2 £F 0 . 



(A .14) 




(A .15) 
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OPERATIONAL AMPLIFIERS AND ANALOG CIRCUITS 



An operational amplifier is a high gain amplifier with 
certain input and output impedances as shown in Figure b.L - 
Analysis of the circuit shows thats 






Z< 



1 + tKz7 + 



where Z f = feedback impedance 
Z^ = input impedance 
A = open loop gain* 



(B.l) 



Note that if A >>(•=—+ 1), this reduces to: 




Since Z f and Z^ are complex quantities, various combinations 
of resistances and capacitances may be used to perform mathe- 
matical operations . Table Bol illustrates the circuits 
which are used to perform the linear operations in the 
analyzing equipment . 




Figo Bol Operational Amplifier Circuit 
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Table B.l Circuitry for Basic Linear Operations Employed in Noise Spectral Analyzer 



APPENDIX C 



ESTIMATION OP k eff AND a 

Estimates of k^ and a can be made by using the 
THERMOS computer code.^ 1 ^ This integral-transport theory 
code, based on the Wigner-Seitz unit cell model, permits 
the calculation ofs 

f^ the thermal utilization 

>7 the average number of neutrons produced 

per thermal neutron absorbed in the fuel 

2 the average macroscopic absorption cross- 
section in the cell 

Combining these values with reasonable estimates 
of s 



7J = age to thermal 
£ = fast fission factor 
p = resonance escape probability 
D = diffusion constant 

B = geometric buckling, 

s 

k , can be estimated bys 
oo 



CD 




(Col) 



where L 2 = — = (diffusion length) 2 
f a 

The lattice investigated in the present studies had 
the following specifications s 
Pitch 2.25 in. 

0.947# U 235 

18.90 gm/cc 



Enrichment 
Fuel Density 
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DgO Purity 99 «4 MPC 

The unit cell dimensions (outer radius minus inner radius) 
are s 

Fuel Oo49l49 cm 

Air gap 0,00762 cm 

Cladding 0.08636 cm 

Moderator 2.4l491 cm 

Using this data, and the procedures outlined in reference 
(14), the following estimates were derived from the THERMOS 
output t 



S = 0.007583 cm” 1 

3 

f = 0.969 
Y( = 1.447 

For this lattice, good estimates for the other quanti» 
ties^^ ares 



D - 8 .75 cm 

p - 0.92 

€ - 1.01 

l * 115,4 cm 2 

B_ = 0.0031 cm” 2 
S 

These values result ins 
k Q0 - 1.30 

k eff - °' 6 5 
JL ~ 4.35 x 10“^ sec 

a y 1 c ^805 sec” 1 

Consulting Figure 1.2 one sees bn. 9 in this lattice, 
required efficiencies for successful noise analysis experi*= 



raents would bes 

for autocorrelations 30# 

for cross-correlation 1.5# 



APPENDIX D 
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COMMENTS ON OPERATION OP EQUIPMENT 
Dol Equipment' 

The equipment assembled to perform the auto spectral 
analysis consists of the following (see Figs. D.l and D„2 ) 0 

1 George A. Philbrick (GAP) MU/DV Duplex Multiplier 
Divider ~ + 100 V input-output* Ser. 0297 

2 GAP Model HK Operational Manifolds* Ser. 069 and 068 

8 GAP K2-XA Operational Amplifiers 

8 GAP K2-P Stabilizing Amplifiers 

8 GAP 2-Unit K-l Modular Assembly Units 

1 Battery Powered AoC. Pre-amplifier utilizing Analog 
Devices Model 106 Operational Amplifier 

2 Associated Specialties Model 2 Power Supply* Ser„ 

748 and 749 

1 Westinghouse Model 6377 Ionization Chamber 
10 No. 467 67 °5 Volt Batteries 

1 Waveforms Model 403B Sine-Square Wave Generator* 

Ser » 7997 

1 Esterline Angus Graphic D.C. Ammeter Model AW* 

Ser o 85339 

Each linear operation is performed using one K2-Xa and 
one K2-P* which are plugged side-by-side into the top of the 
manifold o Each manifold can accommodate 5 pairs of these 
amplifiers . The K-l Modular Assemblies are labeled to desig- 
nate which operation is performed o These assemblies are 
Plugged into the face of the HK manifold so as to complete 
the circuitry for each chopper-stabilised operational ampli- 
fier. Each K-l Assembly is labeled as to input* output* and 
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Preamplifier 

Detector 
Power Supply 




•Multiplier 



K2 



Manifold 
K2-P and K2-XA Pairs 
and K1 Modules 



— -Waveforms 403B 
Oscillator 



sterline 
ec order 



Angus 



Fig. D.l 



Auto- Spectral Density Analysis Equipment 
(Connecting Cables not Shown) 
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Preamplifier 

Detector 
Power Supply 




Multiplier 






K Manifold 
K2-P and K2-XA 
and K1 Modules 



— Waveforms 403B 
Oscillator 




Pairs 



Pig. D. 1 Auto-Spectral Density Analysis Equipment 
(Connecting Cables not Shown) 
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Figo D„ 2 Photograph Showing Filter Module (Left) - 
Preamplifier (Eight ) 
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ground connections . The eight operations performed are, in 
order (see Fig. 4„l)s 

A. Gain of 20 amplification with D.C. blocked 
Bo Gain of 10 or 20 amplification for noise signal 

C. Variable gain (up to 20) amplification for oscillator 
signal 

Do Band pass filter which may be switched to provide either 
narrow bandfo.S cps) or wide band (~10 cps) operation 
Eo Gain of XI, X5, X10, X5©, or X100 swltchable amplifier 
F. Second filter similar to D 
Go Averaging circuit 

The time constant of the averaging circuit may be altered 
by using banana plugs to insert various combinations of 
parallel RC components into the jacks marked FEEDBACK ° 

In order to perform the cross power spectral measurements, 
the following equipment was duplicated 
detector 
preamplifier 

detector battery power supply 

operational amplifiers and modular units to perform 
operations A, B, D, E, and F 
A + 2 5 V GAP MU/DV Multiplier was borrowed from the MIT 
Mechanical Engineering Department to serve as the correlating 
multiplier., A GAP Model R-3QQ Power Supply was also borrowed 
to meet the Increased power requirements. 

D. 2 Equipment Operation 

In order to insure proper operate n, the following had 



to be observed 
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a) Keep equipment well ventilated (e»g» by fan) 

b) Maintain largest signal possible throughout with- 

out exceeding rated signal input limitations 

c) Avoid ground loops 

d) Minimize d.e. amplification at preamp by varying 

resistances in blocking circuit 

e) Avoid preamp oscillations by keeping input cables 

to a minimum length 

These comments supplement the Individual manufacturer^ 
manuals* which have been Included in a booklet Noise Spectrum 
Analyzer* compiled by the author <> 



APPENDIX E 
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PROBABLE ERROR ESTIMATION 

Eol Predicted Error In Amplitude 

The spectral measurements are made using an averaging cir- 
cuit with a time constant = R f C f where R f and are the 
feedback resistance and capacitance in the averaging cir- 
cuit . Using formulas developed by van der Ziel^ 1 ^ and 
Bendat^ 1 ^, the relative standard deviation in cross-power 
spectral density measurements isg 



cr 

_J2L = 

S 



"2 + Nj/S + Ng/S + (Ni/SjCNg/S) 
ttM t> 



1/2 



(Eol) 



where S = strength of signal of interest , 

N 1 2 = stren Sth of undesired signals in channels 1 and 2, 

Af = filter bandwidth, in eps, and 
= averaging circuit time constant 0 
If detector efficiencies are great enough to detect a 
level of correlated noise, S and N would refer to the corre- 
lated and uncorrelated noise, respectively 0 

For most measurements in this work, detector efficiencies 
were so low as to permit measurement of only the uncorrelated 
signal level o For this cases 




Preliminary hohlraum measurements were made with 

Af - Oo88 sec" 1 
X = 100 sec 
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Final hohlraum and lattice auto-correlation measure- 
ments were made with Af =10 sec” 1 

'C = 100 sec 

These result in predicted amplitude errors of 11 *9$ 
and 2.5$, respectively, and are considered to be suffic- 
iently small for this work. If in the future , greater 
accuracy is considered desirable , it may be attained by 
altering the feedback impedance in the averaging circuit 
so as to Increase X = ° 

E .2 Fredicted Error in Frequency 

Error in frequency measurements arise from the finite 
bandwidth of the filters and the lack of accuracy in the 
oscillator dial readings. Frequency resolution due to 
finite bandwidth is + Af , and all equipment dial readings 
are rated to within +2$. 

For example, using a filter with a bandwidth of 10 cps 
to measure the spectral density near a frequency of 300 cps 
would result in an uncertainty in frequency of +LQ cps due 
to finite bandwidth and +6 cps due to dial error. This 
results in a total uncertainty in frequency of +12 cps, or 
about 3.7^o This figure could be reduced by decreasing 
the filter bandwidth. As can be seen, however, changes in 
Af affect both frequency and amplitude uncertainty, thus a 
choice of Af should be governed by desired accuracy in both. 
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